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ABSTRACT 

We analysed Rossi X-ray Timing Explorer observations of the accretion-powered 401 Hz pulsar 
SAX J1808.4— 3658, in order to precisely determine the source distance. While the fluences for the 
five transient outbursts observed from 1996 were constant to within the uncertainties, the outburst 
interval varied signficantly, so that the time-averaged flux (and accretion rate) decreased by around 
40%. By equating the time-averaged X-ray flux with the expected mass transfer rate from gravitational 
radiation, we derived a lower limit on the distance of 3.4 kpc. Combined with an upper limit from 
assuming that the four radius-expansion thermonuclear bursts observed during the 2002 October 
outburst reached at most the Eddington limit for a pure He atmosphere, we found that the probable 
distance range for the source is 3.4-3.6 kpc. The implied inclination, based on the optical/IR properties 
of the counterpart, is z < 30°. 

We compared the properties of the bursts with an ignition model. The time between bursts was 
long enough for hot CNO burning to significantly deplete the accreted hydrogen, so that ignition 
occurred in a pure helium layer underlying a stable hydrogen burning shell. This is the first time 
that this burning regime has been securely observationally identified. The observed energetics of the 
bursts give a mean hydrogen fraction at ignition of {X) w 0.1, and require that the accreted hydrogen 
fraction Xq and the CNO metallicity Zqno are related by Zqno ~ 0.03 {Xq/O.I)'^. We show that 
in this burning regime, a measurement of the burst recurrence time and energetics allows the local 
accretion rate onto the star to be determined independently of the accreted composition, giving a new 
method for estimating the source distance which is in good agreement with our other estimates. 
Subject headings: stars: neutron — X-rays: bursts — pulsars: individual (SAX J1808.4— 3658) — 
stars: distances 



1. INTRODUCTION 

The population of accretion-powered millisecond pul- 
sars (MSPs) has grown rapidly over the last three years, 
with the discovery of six new examples to bring the to- 
tal sample to seven. SAX J1808.4— 3658 remains the 
best-studied of the seven sources, mainly due to the 
fact that it exhibits the most frequent outbursts, every 
2.5 years or so. The source w as first discovered as a 
faint transient with BeppoSAX l|in 't Zand et alllToMl . 
Orbital phase-connection across different outbursts has 
given the most precise determ ination of the orbital pa- 
rameters (Papitto ct al. I2005f) : the mass donor, which 
orbits the neutron star every 2.01 hours, is probably 
a brown dwar f heated by X-ray emis sion during qui- 
escence (Bilds ten &: Chakrabartvll2001|) . The review of 
lwihl^,nds (200^ includes details of the observations un- 
dertaken during each outbur st . Optical obser vations of 
the counterpart V4580 Sgr l|Giles et al.lll999D have led 
to limits on the inclination of cos? = 0.65lg'3?(90% 
confid ence), assuming a distance of 2.5 kpc ( Wan g et al.l 
I2 QQ1[). Th e most recent outburst was in 2005 June 
l|Markwardt et al. 2005; Wijnands ct al. 2005). The neu- 
tron star has been detected in quiescenc e by BeppoSAX 
llSteJIa et all 1200(1 and XMM-Newton l|CamDana et all 
l2002() . and a Chandra observation in outburst revealed 
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the possible presence of a narrow O vii absorption line 
in additi on to neutral absorptio n features attributable to 
the ISM ijCallowav et a,].ll2005^ . 

In this paper, we present a comprehensive study of the 
properties of thermonuclear (type I) X-ray bursts and 
the persistent X-ray emission from the accretion-powered 
MSP SAX J1808.4-3658. We show that the sequence of 
X-ray bursts observed during the 2002 October outburst 
can be understood as a series of thermonuclear flashes 
driven by unstable helium burning in a layer of pure he- 
lium underlying a stable hydrogen burning shell. This 
regime of burning is predicted by X-ray burst models 
when the time between bursts is longer than ~ 1 day, 
sufficient for the hot CNO cycle to exhaust the accreted 
hydrogen and leading to ignition in a pure helium en- 
vironment (Fui lmoto et a l. 1981; Bildsten 1998i) . As far 
as we are aware, however, this regime has not been se- 
curely identified before. Short duration (< 10 s) X-ray 
bursts are observed from Atoll sources at X-ray lumi- 
nosities Lx ^ lO'^^ erg s""'^, and are often referred to as 
"helium-rich" because they do not show the long tails 
characteristic of hydrogen burning. However, these lu- 
minosities are several times larger than those where pure 
helium ignition bursts are expected, and in addition, the 
bursts are not well-described by the standard theory, for 
example having irregular recurrence times and energetics 
that indicate additional burning of fuel between bursts 
fe.g. lBildstellT99ai . 

One of the aims of this study is to compare dif- 
ferent methods used to constrain the distance to 
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TABLE 1 

Outburst properties and calculated distance limits for 
SAX J1808.4-3658 



Start Interval Distance 
Outburst (MJD) (yr)'' Fluence'^ (Fx)" limit (kpc) 



Sep 1996 


50333 


> 0.67 


7.7 ±0.6 


< 36 


1.4 


Apr 1998 


50911 


1.58 


5.2 ±0.5 


10 


2.5 


Jan 2000 


Si 51547 


1.74 


5.4± 1.7 


9.6 


2.6 


Oct 2002 


52559 


2.8 


6.2 ±0.4 


7.1 


3.1 


June 2005 


53522 


2.6 


4.9 ±0.6^* 


5.9 


3.4 



mate is almost independent of the composition of the ac- 
creted material, the main input to the X-ray burst mod- 
els. The plan of the paper is as follows. We present the 
observations and data analysis in §2. In H3.1l we calculate 
the lower limits to the distance obtained by equating the 
long-term time-averaged X-ray flux with the predicted 
mass transfer rate driven by gravitational radiation, ad- 
dressing the issue of long term variations in the measured 
time-averaged persistent flux for the five outbursts ob- 
served so far. In H3.2I we compare the burst properties 
with theoretical ignition models. We first show that the 
burst energetics require a specific correlation between the 
hydrogen mass fraction Xq and the mass fraction of GNO 
elements Zcno in the accreted material. We then use 
this result to show that the measured recurrence times 
and a values imply an accretion rate onto the star that is 
independent of Xq or Zcno , and theref ore g ives a robust 
distance estimate. Finally, in section ti3.3l we estimate 
the distance by equating the peak burst fluxes and the ex- 
pected Eddington fluxes reached during radius-expansion 
bursts. We conclude in §4. 

2. OBSERVATIONS AND ANALYSIS 

We analysed archival observations of 
SAX J1808.4-3658 made with the Rossi X-ray 
Timing Exp lorer (RXTE). Th e proportional counter 
array (PGA; iJahoda et allllOQfiH aboard RXTE consists 
of five identical proportional counter units (PGUs) 
sensitive to photons with energies between 2-60 keV, 
with a field of view of 1° and a total collecting area 
of « 6000 cm^. In order to characterise the pulsar 
spectra up to w 100 keV, we also used data from 
the High-Energy X-ra y Timing Experiment (HEXTE; 
Rothschild ct alT I1998D . which consists of two clusters 
of four NaI(Tl)/GsI(Na) phoswich scintillation counters 
also with a field of view of 1°, sensitive to photons with 
energies between 15-250 keV and with a total collecting 
area of 1600 cm^. 

We analysed the data with lheasoft version 5.3, re- 
leased 2003 November 17. We extracted PGA Standard-2 
data over each observation (excluding any bursts) , sepa- 
rately for each PGU. We fitted absorbed blackbody plus 
powerlaw models in the energy range 2.5-25 keV, inte- 
grating the model to derive the persistent source flux Fp. 
For selected (preferably long) observations we made ab- 
sorbed Gomptonization ("comptt" in XSPEC) model fits 
to combined PGA and HEXTE data, plus a Gaussian 
describing Fe Ka emission around 6.4 keV, where neces- 
sary. We then generated an idealised response with 200 
logarithmically-spaced energy bins between 0.1-200 keV 
and integrated the model flux over this range, using 
the ratio of the 0.1-200 keV to 2.5-25 keV flux as a 
correction Cboi in estimating the bolometric flux from 
each observation of each source. The X-ray colors for 
SAX J1808.4— 3658 were relatively constant throughout 
each outburst, and so we adopted a constant value of 
Cboi = 2.12. We used measurements of the burst prop- 
erties tabled in the catalog of Galloway et al. (2006, in 
preparation; hereafter G06). These include the start time 
(measured as the time the burst flux first exceeds 25% 
of the maximum), bolometric fluence Eh, and the ratio 
of the integrated persistent flux to the fluence a. The 
latter parameter is subject to systematic errors due to 
intermediate bursts which may be missed during source 



''The epoch for the outburst prior to the first known is assumed 
to be earlier than the first ASM measurements (1996 January 6, 
or MJD 50088). 

''Estimated bolometric fluence in units of 10~^ ergscm"^ 

'^Estimated mean bolometric flux in units of 
10~^^ ergscm~^s~'. 

"^The fluence for the June 2005 outburst was estimated primarily 
from the ASM observations, since no public PCA data from the 
first 20 d of the outburst were available 



SAX J1808.4-3658 and other MSPs. Two principal 
methods have been used to date. First, three of the seven 
sources have shown thermonuclear X-ray bursts, which 
have a maximum peak flux believed t o correspond to the 
Eddington limit at t he surface (e.s. ILewin et alJ 119931: 
IKuulkers et al.ll2'003() . The uncertainty in distances de- 
termined by this method arises from uncertainty in the 
mean hydrogen fraction at the photosphere, and the neu- 
tron star mass. The distance to SAX J1808.4— 3658 was 
initally estimated at around 4 kpc based on the pe ak flux 
of a thermonuclear burst (jin 't Zand et al.lll998j) . The 
likely distance range was later revised to between 2.5 and 
3.3 kpc; the decrease from 4 kpc resulted from a more 
accurate analysis of an additional burst not detected in 
the original data, while the range arises from the possible 
variation of the Eddi ngton limit betweeen a n H-rich and 
pure He atmosphere l|in 't Zand et aljr200lD . 

Second, the distance to a transient binary can be 
constrained by equating the time-averaged X-ray flux 
(estimated from the outburst fluence divided by the 
time since the previous outburst) and the flux predicted 
from the expected mass transfer rate due to gravita- 
tional radiation Mgr., which depends on the orbital pe- 
riod and the response of the mass donor to the loss 
of material (the mass-radius re lation; see for example 
iVerbunt fc van den Heuveilll995D . Because pulse timing 
gives the projected semimajor axis ax sini [i is the incli- 
nation of the orbit) , only a lower limit on the companion 
mass is possible, which gives a lower limit to the accre- 
tion rate and the distance. An important issue is the long 
term X-ray flux behavior of the MSPs. The mass trans- 
fer rates in many non-pulsing low-mass X-ray binaries 
(LMXBs) vary on time scales comparable to MSP out- 
burst intervals, and it is possible that the same is true for 
the MSPs. Presently this remains an open question be- 
cause the MSPs have uniformly low accretion rates, and 
are transient, so that for most sources only one outburst 
has been observed. 

In this paper we also show that models of the sequence 
of X-ray bursts from SAX J1808.4-3658 give a third es- 
timate of the accretion rate onto the star and therefore 
distance to the source. Remarkably, this distance esti- 
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Fig. 1. — Predicted burst recurrence times according to the 
ignition model for solar metallicity {Zqt^q = 0.016) in the accreted 
material. The flux from the crust Qj, is fixed at 0.3 MeV/nucleon, 
and we show curves for a range of hydrogen fraction Xq in the 
accreted material. The s quare symbols show th e prediction for the 



time-dependent model of lWooslev et aLI 1200' 



occultations occuring every 90 min orbit. 

We compared the observed pr operties of the bursts 
with the ignition calculations of lCumming~ Bildstcn 
l)200Cl|) . to which we refer the reader for further de- 
tails. We assumed a 1.4 M© neutron star with ra- 
dius i?NS — 10 giving a surface gravity g = 
{GMj^s/Rm)i^ + ^) 2.45 x 10" cm s'^ where z is 
the redshift at the NS surface. For our choice of mass 
and radius, 1 + z = (1 - 2GMj^s/ Rnsc^)~'^/^ = 1.31 
(cf. iCottam et all ' 20021. For a given accretion rate m 
(parameterised as a fraction of the Eddington accretion 
rate mEdd = 8.8 x 10^'' gcm"^) we calculated the tem- 
perature profile of the accumulating layer of hydrogen 
and helium, adjusting its thickness until the condition 
for a thermal runaway is met at the base. The tempera- 
ture is mostly set by hydrogen burning via the hot CNO 
cycle, and therefore the CNO mass fraction Zqno, which 
we refer to as the metallicity. The model also includes 
a flux from the crust Qh. To calculate the burst energy, 
we assume complete burning of the H/He fuel layer, and 
that the accreted material covers the whole surface of the 
star. The total energy is then 47ri?^y(5nuc/(l + z), where 
y is the ignition column depth, and Qnuc is the energy per 
gram from nuclear burning. We write Qnuc = 1.6-1-4 {X) 
MeV per nucleon, where {X) is the mass-weighted mean 
hydrogen fraction at ignition. This expression for Qnuc 
assumes « 35 % energy loss due to n eutrinos during the 
rp process (e.g. lFuiimoto et aPHQSTl) . and gives 4.4 MeV 
per nucleon for approximately solar hydrogen abundance 
{X = 0.7). 

Figure n] shows the predicted burst rate as a function 
of m for different accreted hydrogen mass fractions. The 
open squares are the results of the time-dependent cal- 



culations of burst sequences by iWooslev et al.l l)2004(l . 
At low m < 0.05 TTiEdd, the agreement between the 
two calculations is good. At higher accretion rates, 
and therefore high burst rates, the residual thermal en- 
ergy from previous bursts becomes an important fac- 
tor in determining the i gnition point, giving ri se to 
the difference betw e en thelCumming fc B ildstcn F2000f) 
and IWooslev et alJ l|2004|) calcul ations (see the discus- 
sion and comparison in §5 of IWooslev et all |2004). 
Fortunately, the range of accretion rates spanned by 
SAX J1808.4-3658 is almost always < 0.05 mEdd- 

3. RESULTS 

3.1. The long-term mean accretion rate and 
corresponding distance lower limit 

We analysed all the available outbursts from 
SAX J1808.4— 3658 to calculate the time-averaged flux 
from interval to interval. The timc-averaged accretion 
rate driven by angular moment um loss by gravitational 
radia tion from the binary is ijBildsten & Chakrabartvl 
12001 
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1.4Mo/ I o I,. / (1) 

where Mc is the minimum companion mass, and 
Porb the binary orbital period. Assuming conserva- 
tive mass transfer, we equated the measured time- 
averaged bolometric X-ray flux (Fx) and the predicted 
flux {GMms / Rns)Mq's./ 'ind^ to derive a lower hmit on 
the distance of 

, '1/2 

d > 2.6 kpc ' 
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10 1° ergs cm ^ s ^ 
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-4/3 



NS 



1.4 Mp 







5/6 
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where we choose the minimum companion mass 
Mc ~ O.O43M0, appropriate for the mass func- 
tion of SAX J1808.4-3658 of 3.78 x IQ-^A/q 
l|Chakrabartv fc Morgai3ll998j) and assuming a neutron 
star mass of 1.4A^0. 

SAX J1808.4— 3658 is the only accretion-powered 
MSP in which multiple outbursts have been observed 
with RXTE/PCA, although outbursts have been de- 
tected retroactively in the long-term ASM flux histories 
of both XTE J 17 51-30 5 and ICR 100 291-1-5934 
fMarkwar dTeTan l200l IGallowav etUI |2005). 
We estimated the fluence for each outburst of 
SAX J1808.4-3658 using public RXTE PCA and 
ASM measurements, shown in Figure |2| Generally, the 
PCA observations (which provide the best measure of 
the source flux) are triggered after the outburst has 
commenced, and often do not cover the entire outburst, 
unlike the ASM measurements. For this reason we 
integrated the ASM intensities over the parts of the 
outbursts not covered by the PCA measurements, and 
determined a cross-calibration between the PCA flux 
and the ASM intensity, as follows, for each PCA flux 
measurement we selected the closest 1-day average of 
the ASM 2-10 keV intensities, and performed a linear 
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TABLE 2 

Properties of X-ray bursts from SAX J1808.4— 3658 



Burst start time Burst start time Peak flux Fluence 

(UT) (MJD) (10^9 ergscm-^s-i) (IQ-'' ergscm^^) 

1 2002 Oct 15 09:55:37 52562.4136 163 ± 3 2.071 ± 0.015 

2 2002 Oct 17 07:19:24 52564.3051 177 ± 3 2.139 ±0.012 

3 2002 Oct 18 04:25:20 52565.1843 179 ± 3 2.436 ± 0.012 148 ± 3 

4 2002 Oct 19 10:14:33 52566.4268 177 ± 3 2.716 ±0.015 167 ± 3 



''The a values are calculated between pairs of bursts, excluding the first two since we expect that 
intermediate bursts may have been missed. 
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Fig. 2.— RXTE ASM and PGA X-ray flux measurements 
throughout the five outbursts of SAX J1808.4— 3658, from top to 
bottom: 1996 September, 1998 April, 2000 January, 2002 October 
and 2005 June. The gray boxes show the Icr confidence regions 
for the 1-d averaged ASM intensity estimates (left-hand y-axis). 
The open squares show the PCA 2.5-25 keV fiux measurements 
(right-hand j/-axis), with error bars indicating the Itr uncertainties 
(typically smaller than the symbols). The times for each outburst 
have been shifted relative to the estimated outburst start times, 
listed in Table [Tl Note that additional proprietary PCA observa- 
tions (PI: Wijnands) were made in the first 20 days of the 2005 
June outburst. 



fit. The fluences derived by this method, scaled to give 
the estimated bolometric values, are listed in Table 
We propagated the errors from the errors on individual 
ASM/PGA measurements. We note that these values 
are roughly consistent with prior estimates, e.g. for the 
1998 April outburst iGilfanov etaP l)1998(l estimated 



4.2 X 10 ^ ergs cm ^ (3-150 keV), whereas we estimate 
(5.2 ±0.5) X 10-3 ergscm-2 (extrapolated 0.1-200 keV). 

The derived fluences were constant to within the errors, 
with variation only at the 2.6(t level. However, the inter- 
val between outbursts increased after each of the the first 
three, from 1.58 yr between the 1996 September and 1998 
April outbursts, to 2.8 yr between the 2000 January and 
2002 October outbursts. Although the interval to the 
2005 June outburst was slightly less than the previous in- 
terval, the integrated fluence was the lowest of all the out- 
bursts. As a result, the time-averaged X-ray flux dropped 
by 40% between 1996 and 2002. Since the derived dis- 
tance limit depends on the mean flux only to the —1/2 
power (equation [21 , the derived distance limits varied 
only by 30%, up to a maximum of 3.4 kpc (Table^. The 
mean X-ray flux, based on the four outbursts observed 
since 1996 September, was 8x 10"^^ ergs cm~^ s"^, which 
(for a distance of 3.6 kpc) corresponds to approximately 
10-" Moyr-^ 

3.2. Comparison of burst properties to ignition models 

We now compare the properties of thermonuclear 
bursts from SAX J1808.4— 3658 with predictions from 
theoretical ignition models. Observed properties of the 
bursts are listed in Table 2. The bursts were quite homo- 
geneous, with standard devation of the fluences of only 
about 12%. The bursts were all observed within a 100 hr 
interval just after the outburst peak, and the last three 
were separated by 21.1 and 29.8 hr, respectively. This 
increase in burst interval was observed as the 2.5-25 keV 
flux decreased from 2.3 to 1.9 x 10"^ ergscm"^ s"^. Gaps 
between the observations with bursts make it possible 
that intermediate events were missed. Assuming that 
was not the case, the estimated a-values for the second 
two intervals were a — 148±3, 167±3 respectively. These 
values, along with the fast rise times (« 0.5 s) suggest 
pure He fuel. This is supported by the lightcurves of the 
X-ray bursts, an example of which is shown in Figure 
The long « 10 second period of constant (presumably 
Eddington) luminosity followed by a deca y over » 10 sec- 
onds is remarkably similar to the results of lWooslev et alJ 
l|2004'l for bursts which ignite helium in a pure helium 
environment (compare with their Fig. 26). As we will 
show, the recurrence times and energetics of the bursts 
also suggest that helium ignition is occurring in a pure 
helium layer beneath steady hydrogen burning shell. 

3.2.1. Pure-helium triggered X-ray bursts 

Before presenting our detailed fits to the ignition mod- 
els, we first make some analytic arguments which help 
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Fig. 3. — The fourth, and most intense, burst observed from 
SAX J1808.4-3658 by RXTE during the 2002 October outburst. 
Top panel The estimated bolometric flux is plotted as a histogram, 
with error bars indicating the la uncertainties; at the burst peak, 
only short intervals of data were obtained due to telemetry buffer 
overruns. The distance is assumed to be 3.6 kpc. Middle panel The 
blackbody temperature fcT^b- Bottom panel The blackbody radius 
Rhh- Note the anticorrelated behavior of fcT^jj and ij^b during the 
burst peak. The circle marks the time when the peak burst flux 
was achieved. 



to understand the observed properties of the bursts and 
the results of our fits. First, consider the burst ener- 
getics. The observed values of a « 150 imply that the 
nuclear energy release per gram is Qnuc = c^z/a « 2 
MeV per nucleon, where we take the gravitational red- 
shift 1 + z = 1.31. Using the relation Qnuc = 1-6 -I- 4(X) 
MeV per nucleon gives an average hydrogen mass frac- 
tion {X) « 0.1 in the burning layer immediately prior to 
ignition. 

The low mean hydrogen fraction at ignition comes 
about because hot CNO hydrogen burning leads to 
steady depletion of the hydrogen as the H/He layer 
accumulates. The time to burn all the hydrogen in 
a given fluid element depends on the initial hydrogen 
fraction Xo and the mass fraction of CNO elements 
•^CNO- In the observer's frame, this time is ^burn = 
38 hrs Xo (Zcno/0.01)-1(1 + z)/1.31. The bursts from 
SAX J1808.4— 3658 have recurrence times At « 1 day 
and {X) « 0.1. This immediately constrains .^cno 
and Xq: only particular combinations of these two pa- 
rameters will lead to an average hydrogen mass frac- 
tion w 0.1 for the observed recurrence time of « 1 day. 
The average hydrogen fraction in the layer at ignition is 
(X) = (Xo/2)(tburn/At), where we have used the fact 
that the hydrogen fraction d ecreases linearly with depth 
(|Cumming fc Bildstenll2000|) . Setting {X) and At to the 
observed values, we find 



Xo = 0.36 



■^CNO 

0.01 



1/2 



At 
1 day 



1/2 



\0.1 



1/2 



l + z 
1.31 



-1/2 



(3) 



is the required relation between Xq and .Z^cno to repro- 
duce the measured values of At and {X)^. 

To go further, we need to understand the ignition col- 
umn depth. An example of a temperature profile in the 
layer for At > tburn is given in Figure 1(b) of Gum- 
ming & Bildsten (2000). The temperature rises through 
the hydrogen burning layer as the hydrogen burns, but 
the layer of pure helium which accumulates beneath it 
is close to isothermal. Therefore the temperature at the 
base of the hydrogen burning layer is a good approxima- 
tion of the ignition temperature. The temperature pro- 
file in the hydrogen burning layer is given by dT/dy — 
iFn/AacT^, where k is the opacity, and y is the column 
depth (units of g cm~^). The flux is F = eniyd — y) 
where yj^ = uiXqEh / is the depletion depth for hy- 
drogen. The energy release in the hot CNO cycle is 
Eh = 6.0 X 10^* erg g~^; the heating rate from the hot 
CNO cycle is = 5.8 x lO^^ erg g-^ s~i (Zcno/0.01). 
Integrating from the surface, we find the temperature at 
the base of the hydrogen layer is given by^ 

4 _ iK ra'XlEl Xl 

Jfc — °^ • \^) 

lac en ^cno 

But we have already seen that the energetics fix the ratio 
Xg/ZcNO (eq. 0), and so the temperature Tj, is inde- 
pendent of Xq or Z'cnOj 

1/4 

/ .^K TT f I I ( /\ )/\/ \ 

J?) — 
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\ ac(l -I- z) 



:2.1 X lO** K 



1/2 



1/4 



/ At 
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({XT 

V 0.1 J V24 h, 
as a typical value for 



^O.lTOEdd, 

where we set k — 0.05 cm^ 
opacity and z — 0.31. 

Next we assume that the helium layer is isothermal, 
so that Tb is also the temperature at the helium ignition 
depth. It can be shown that this is a good assumption 
as long as the flux heating the helium layer from below 
is Qb < Eh{X) « 0.6 MeV per nucleon for {X) of 0.1, 
which is satisfied for the expected values of Qb ~ 0.1 
MeV per nucleon (e.g. Brown 2000). A good fit to the 
helium ignition curve of ICumming fc BildstenI l|2000() is 
yb = lOi" g cm-2 (Tb/1.1 x 10^ K)-^ (for > 2 x 
10* g cm^^). Substituting equation Q for Tb into this 



fit, and solving for m 
m = 0.057 TTiEdd 



yb{l + z)/At gives 
At \-''/'^ 

I 0.1 



24 h 



(6) 



^ Although not shown explicitly, the normalization of equation 
13 is quite sensitive to the assumed redshift through the factor < 
X >, since for a given a, Qnuc oc z. For example, a larger choice 
for the redshift of z = 0.42, which gives (X) = 0.25, shifts the 
relation in equationl^so that it admits both solar values for Xq and 
•^CNO- However, as we describe in §3.2.2, our detailed flts, which 
account for the variations in burst properties from burst to burst, 
constrain the relation between Xq and .Zcno to be approximately 
that given by equation^] or in other words imply that z 0.3. The 
resulting constraints on mass and radius are described in §3.2.2 (see 
discussion below eq. 13). 

^ The base temperature Tj, has some dependence on composi- 
tion via the opacity. However, this is a small effect. Varying the 
hydrogen fraction or metallicity from zero to solar values, we flnd 
that the resulting changes in opacity are < 50%, and the changes 
in temperature (eq. |^ and inferred accretion rate (eq. |S} smaller 
still because of the weak dependences Tf, oc k^^"^ and rh oc k"^^^^ . 
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Fig. 4. — Comparison of model results and observed burst 
properties for SAX J1808.4-3658, for Qb = 0.325 MeV/nucleon, 
ZcNO = 0.02, Xo = 0.54 (x^ = 6.2 for 6 dof). The top panel 
shows the inferred bolometric flux [open squares) and the inter- 
polation used to estimate the integrated m between the bursts 
{dashed line). The open triangles are the observed burst times, 
and the open diamonds indicate the error (observed-predicted) in 
the model burst times (relative to (2). Lower left The observed 
and predicted fiuence, the latter rescaled by r2; the dashed line 
indicates a 1:1 correspondence. Lower right The observed and pre- 
dicted a- values, the latter rescaled by r:j. The measurement for 
burst 2 is a mean assuming that two bursts were missed between 
burst 1 and 2. 



where we again set k = 0.05 cm^ g"-'^. 

To summarize, we have found that given a sequence 
of pure helium ignition bursts with measured At and 
{X) (or equivalently At and a), Xq and .^cno rnust 
be related according to equation |(2Jl, and the accretion 
rate is determined by equation ^ independent of Xq 
or ZcNO- Writing Fx = mQgrav(-R/(i)^/(l + z) where 
Qgrav — c^z/ll + z) GM / R is the gravitational energy 
release per gram, and solving for the distance, we find 



d = 4.15 i?6 
At 



F 



X 



5 X 10 ^ erg cm ^ s 



24 h 



-13/44 



V 0.1 



-9/44 



kpc 



(7) 



which is independent of Xq and Zq^q . 



3.2.2. Fitting procedure and results 

We can improve on the estimates of the previous sec- 
tion by firstly fitting the observations to the full ignition 
models (which, for example, include accurate opacity val- 
ues as a function of depth in the layer), and also by fit- 



ting for the change of burst recurrence time with chang- 
ing accretion rate during the outburst. We generated 
sequences of bursts based on the inferred rh history of 
SAX J1808.4-3658 during the 2002 outburst as follows. 
First, for a given set of model parameters (Qf,, ^cno and 
Xq) we identified the values of m for which the recurrence 
time was closest to that of the third and fourth bursts (as- 
suming no intermediate bursts occurred), At2 = 21.10 hr 
and Ata — 29.82 hr. From these m2, m-i values we cal- 
culated the ratio of the persistent fiuxes to the inferred 



m, I.e. 



ri 



(8) 



/ Thj 
\ Cho\Fp, 

where Fp i is the mean 2.5-25 keV persistent flux (in 
units of 10~^ ergscm~^s~^, from cubic spline interpola- 
tion) between burst (z — 1) and i. The value of ri will 
obviously depend upon the distance to the source, as well 
as the neutron star mass and radius, but at this stage we 
assume no particular values for these parameters, only 
that the relationship between Fp and m is constant over 
the interval where the bursts occurred. 

We converted the interpolated flux history of 
SAX J1808.4— 3658 to an m history by scaling by ri, and 
generated sequences of bursts forward in time from the 
first burst using the ignition model. We then matched 
simulated bursts predicted to occur at approximately 
the same time as the observed bursts, and compared 
their properties. Specifically, we compared the model 
predicted variation of the fiuence and a-values by cal- 
culating two additional ratios r2 = (where 

ttprod = 290/Qnuc.pred IS the predicted value from the 
model, and a — FpC\jo\At / is the observed value), and 



f3 



(where Lb,pred is the predicted burst lu- 
minosity in 10'^^ ergs s~^). We scaled aprod and Lb^prcd by 
the factors r2 and rs, respectively, in order to compare 
with the observed values. Effectively, we try to match 
the observed variation of Ej) and a with m, under the 
assumption that the average offset between the observed 
and predicted values is a function of the system param- 
eters (see §3.2.3). In order to quantify the agreement 
with the observed values, we calculated an overall for 
the simulated burst properties, including predicted burst 
times for the second through fourth bursts; 



^2_^ / a - r2apred j ' 



E 



t, - 1 



z.prcd 



(9) 



where we adopted a representative uncertainty for the 
burst arrival times of Ct = 10 min chosen so that the 
contribution of the last term in equation (jSJ was com- 
parable to the other terms in the sum (the burst start 
times are actually measured to a precision of less than a 
second). 

We carried out the simulation for a range of in- 
put parameters defined by grids of Qb between 0.175- 
0.5 MeV/nucleon (in steps of 0.025), Zcno = 0.006- 
0.034 (steps of 0.002) and Xq = 0.28-0.74 (steps of 0.02), 
and calculated the for each grid point. The burst 
times were reproduced well over a range of the input pa- 
rameters, with an error as small as just 0.2 hr rms. The 
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variation in Ef, and a was less well fitted, giving an over- 
all Xmin — 49. With three burst arrival times, three a 
values and four measured fluences, and three parame- 
ters in the fit {Qb, -^cno and Xq) the number of degrees 
of freedom is 6, so that the best fit is not statistically 
acceptable. 

G06 estimated the errors on the fluence for the bursts 
by propagating the uncertainties in the measured flux in 
each time bin. However, during the 6-8 s PRE episode 
in each burst the X-ray flux was between 1.6 x 10~^ and 
1.8 X 10"^ ergscm~^s~^, so that the onboard teleme- 
try buffer was filled more rapidly than it could be read 
out (once per second; Fig. OJ. As a result, we have in- 
complete information about the flux at the peak of the 
bursts, which makes it likely that the uncertainty on 
the fluences (« 0.6% on average) was underestimated. 
The variation of the flux at the peak of the burst of 
~ 10~® ergscm~^s~^, combined with the typical obser- 
vational duty cycle there (due to the gaps) of « 30% 
suggests that the actual fluence uncertainty is around a 
factor of 2 higher. By increasing the error on the fluences 
by this factor, and propagating the uncertainties for the 
a calculations, we obtained a value which was still sig- 
nificantly larger than the number of degrees of freedom. 
Increasing the errors on the fluence by an additional fac- 
tor of two was sufficient to achieve a minimum of 6.2, 
comparable to the number of degrees of freedom. We 
compare the measured burst properties and the best set 
of model results graphically in Fig. 01 

We also estimated from the grid results the confi- 
dence region for each of the three parameters. The low- 
est values of were consistently achieved for Qt = 
0.3 MeV/nucleon. In order to illustrate graphically the 
combined confidence region for the other two parameters, 
we fiattened the grid along the Qb axis to plot Figure 
m We adopt as a confidence level the probability cor- 
responding to a la deviation, for which the customary 
confidence region has Ax^ = — < 2 30 (f or two 
parameters of interest; e.g. iLampton et al.llT976() . We 
found the best agreement of the measured burst prop- 
erties to the ignition model within a diagonal strip in 
^CNO'-'^o space. This correlation is in good agreement 
with equation which follows from consideration of 
the burst energetics. We could not obtain an accept- 
able model fit for Xq < 0.3. The best solution for solar 
■^CNO = 0.016 required Xq — 0.5; in contrast, for a so- 
lar Xq — 0.7, the corresponding ^cno was substantially 
above the solar value at 0.03. 

As a consistency check we compared the predicted 
times of bursts which were not detected in the RXTE 
observations, with the time intervals in which RXTE was 
observing the source. We predicted two additional bursts 
on October 16 (Fig. 0J, between the first two bursts de- 
tected by RXTE; none of the model runs which agreed 
well with the observations (i.e. Ax^ < 3.51) predicted 
times for these bursts which fell within the observations 
on that day. We also examined the seventh burst in 
the series, predicted to follow the last burst detected by 
RXTE on October 19. The set of models that agreed 
well with the observations predicted times for the seventh 
burst between MJD 52568.398 and 52568.601 (October 
21); just five predicted times fell within the time interval 
for one of the observations on that day, compared to 273 
which fell within the gap. Thus, we estimate the prob- 
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Fig. 5. — Contours of as a function of Zq^q, Xq for the igni- 
tion model fits to the 2002 Octo ber b urst sequence observed from 
SAX J1808. 4-3658, described in JO The dotted Hue corresponds 
to solar metallicity, while the dashed line indicates solar H-fraction. 
The best solution for each grid value of Zcno is indicated by the 
open square. 



ability of having observed this burst as 0.018, which is 
consistent with the lack of detection by RXTE. 

The ignition models also allow us to estimate the mean 
H-fraction at ignition for each of the bursts. Over the la 
confidence region in Zc]<io-Xq space, the model-predicted 
mean hydrogen fraction at ignition {X) for each burst 
varied by ^ 10%. For the four bursts observed by RXTE, 
we find {X) = 0.191 ± 0.016, 0.176 ± 0.015, 0.133 ± 0.011 
and 0.095 ± 0.008, respectively. The decrease in (X) is 
consistent with the steadily increasing burst recurrence 
time, which allows more of the accreted H to be removed 
by steady burning. 

3.2.3. Distance estimate 

Based on our simulation results, we attempted to con- 
strain the distance to SAX J1808.4— 3658. The accretion 
rate (per unit area) is 



^pCbol(l + z) 



Q 



grav 



:6.7 X 10^ 

Mns 



10 ^ ergs cm 

1 



10 kpc 



1.4 M, 







1 



z 



1.31 J VIO km 



gcm~2(5ltf) 



where d is the distance to SAX J1808.4— 3658, and 
Qgrav — c^z/{l + z) « GA/nsZ-^ns IS the energy released 
per nucleon during accretion. Then 



ri= 0.0762 
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Fig. 6. — Peak burst flux as a function of the mean hydrogen 
mass fraction {X) inferred from the ignition model, for the four 
bursts observed by RXTE from SAX J1808.4-3658 during the 
2002 October outburst. Error bars on the flux indicate the la 
uncertainties; for {X), the range spanned by the set of model results 
for which the values are within the Icr confidence limits in Fig. 



1 + z 



Rns 



(11) 

1.31 y VlO km/ ^ ' 

Most of the terms in equation 1111 are expected to be of 
order unity. Thus, we adopt di = lO-yri/0.0762 kpc as 
a rough estimate of the distance; we found a range of 
3.1-3.5 kpc over the set of model parameters Qb, Zcno 
and Xo for which the value was within th e la l imit 
for 3 parameters (Ax^ = 3.51: .Lamoton et al.ill97(f) . 
Similarly, we calculate 



r3 = 82.6 



10 kpc 



-Rns 
10 km 



1.31 



^5nuc,prcd 



(12) 



so that we obtain a second estimate of the distance — 
10-\/82.6/r3 kpc based on the observed and predicted 
fluences. Over the Icr confidence range of ignition model 
parameter space we find a distance range of 3.6-3.8 kpc. 
That the ranges of di and ^3 do not coincide indicates 
that at least one of the terms in equations ^] and [T^ 
deviates significantly from 1. Because the i?NS term is 
the only one common to both expressions, and it appears 
in both distance expressions with exponents of the same 
sign, we infer that no single term can give rise to the 
disagreement. 

The calculated values of r2 can give an indication of 
which parameters may vary from our estimated values. 
We derive 



r2 = 0.668 



M 



1AM, 



R 



10 km 



0„ 



Qnuc,prcd 



(13) 



and we found over the 1 — ct confidence region r2/0.668 
varies between 1.56 and 1.79. This suggests that the 
neutron-star mass must exceed IAMq, or that the radius 
must be < 10 km or that Qnuc < Qnuc.pred (or some 
combination of the three). A larger Mns > IAMq would 
explain the lower value of ri (equation^} leading to the 
distance estimate di < d^; similarly, Qnuc < Qnuc.prcd 
would lead to an inflated value of d^. On the other hand, 
a smaller i?NS < 10 km would drive both di and ^3 to 
even greater disagreement. Thus, we conclude that the 
likely explanation for the disagreement between di and 
ds is that Mns > IAMq and Qnuc < Qnuc.prod, and that 
the distance range is between 3.1-3.8 kpc. 

3.3. Distances from thermonuclear burst peak fluxes 

Presently, the most commonly-used method to deter- 
mine the distance to a bursting LXMB is to measure 
the peak flux of thermonuclear bursts which exhibit pho- 
tospheric radius-expansion. Such bursts exhibit a lo- 
calized maximum in the fitted blackbody radius in the 
first few seconds of the burst, accompanied by a local 
minimum in the blackbody temperature. Several au- 
thors have shown that the peak flux of radius-expansion 
bursts for sources with known dist a nces i s a rough stan- 
dard candle; e.g. iKuulkers et alJ l)2003D found that a 
sample of radius-expansion bursts from globular cluster 
LMXBs almost all reached a common peak luminosity 
of (3.79 ±0.15) X 10"^^ ergss~^. This empi rical value is 
some what higher than predicted by theory ijLewin et 
IT99I 



87rGmpMNSc[l + (axTc 



NO. 861 



(Tt(1 
= 2.7 X 10^* 

^l + z{R) 



X)[l + z{R)] 
Mns \ 1 + (aTTe)0-86 



1.31 



1.4M0 

1 

ergs s 



(1 + ^) 

(14) 



where To is the effective temperature of the atmosphere, 
ax is a coefficient parametrizing the temperature de- 
pendence of the ele ctron scattering opacity (~ 2.2 x 
10^^ K^^; iLewin e t al. 1993) and X is the mass frac- 
tion of hydrogen in the atmosphere at ignition. The 
final factor in square brackets represents the gravita- 
tional redshift at the photosphere 1 + z{R) = (1 — 
2GMns/-Rc2)-i/2. The photosphere may be elevated sig- 
nificantly above the NS surface (i.e. R > i?Ns)i in which 
case this is a small correction. 

All four bursts from SAX J1808.4-3658 observed 
by RXTE during the October 2002 outburst showed 
unambiguous indications of radius expansion (see also 
fChakrabartv ct al. 2003). These bursts all exhibited a 
local maximum in fitted blackbody radius and local min- 
imum in blackbody temperature (Fig. immediately 
prior to the time at which maximum flux was reached. 
The average peak flux was (174±8)xl0~® ergscm~^s~^, 
indicating a distance oi d = 2.8 (3.6) kpc for X — 
0.7 (0.0) in equation [H (see also G06). These val- 
ues are approximately the same as those derived by 
lin 't Zand et al.l 1)200 If) from the bursts observed first by 
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BeppoSAX during the 1996 September outburst. The 
corresponding d istance assuming the empirical value of 
iKuulkers etall l)2003() is 4.3 kpc. The best agreement 
with our distance estimate from §3.1 is for X « at the 
photosphere. Without a detailed model of the expansion 
of the outer layers, it is not possible to say how the com- 
position at the photosphere relates to the mean hydrogen 
fraction in the layer. Nonetheless, we note that there is a 
rough anticorrelation between the peak flux of the burst 
and {X) from our ignition model fits. This is shown in 
Figure El 

4. DISCUSSION 

We have analysed the X-ray flux evolution during five 
outbursts of SAX J1808.4— 3658, as well as the sequence 
of thermonuclear (type I) X-ray bursts observed during 
the 2002 October outburst. The large a- values estimated 
for the bursts indicates low mean H-fraction (X) « 0.1 
at ignition, which analytically implies a common con- 
straint on the accreted H-fraction Xq and metallicity 
^CNO- We confirmed this result by comparing the mea- 
sured burst properties and recurre nce times with predic- 
tions by the burst ignition model of lCumming fc BildstenI 
(pOOO). We found that the inferred burst recurrence 
times were sufficiently long that the helium ignition oc- 
curs in a pure helium environment, the accreted hydro- 
gen having completely burned away between bursts by 
the hot CNO cycle. While the model comparison admits 
a range in both .^cno and Xq, it excludes solar values 
for both. If the metallicity is approximately solar, then 
the hydrogen fraction in the accreted material Xq must 
be significantly sub-solar; on the other hand, if Xq is 
approximately at the solar value, then ^cno must be 
significantly above the solar value. Many stars in the 
Galaxy have metallicities twice solar, and the progenitor 
of SAX J1808.4-3658 may have been a 2-3 M© star. 
The metallicities of other stars in the neighborhood may 
provide circumstantial evidence in support of super-solar 
values for SAX J1808.4-3658, although this would be far 
from conclusive. 

We used our fits to the burst models to constrain the 
distance to the source. Our analytic estimates suggest 
that the local accretion rate onto the star, and therefore 
the distance, is determined independently of the accreted 
composition in this burst regime. Our model fits support 
this, giving a distance range of 3.1-3.8 kpc across the grid 
of the parameters Xq, Zcno^ and Qb- We also used fact 
that the bursts showed photospheric radius expansion to 
constrain the distance, assuming that the peak burst flux 
was equal to the local Eddington luminosity. This gives 
a distance of 3.6 kpc for a pure helium photosphere. 

Finally, we obtained a constraint on the distance based 
on measurements of the fiuence and recurrence time for 
each of the outbursts observed since 1996 September. 
The fluences were approximately constant between out- 
bursts, but the increase in the interval between outbursts 
indicates that the long-term accretion rate is steadily in- 
creasing. We derived a lower limit on the distance based 
on the assumption that the long-term mass transfer rate 
is driven by gravitational radiation in the binary. Since 
the long-term average flux is decreasing, this limit has 
increased with each successive outburst, to a maximum 
of 3.4 kpc for the interval prior to the 2005 June out- 
burst. This limit excludes all but the largest values in 



the range derived from the peak burst flux. By com- 
bining this limit and the estimated range based on the 
peak burst flux, we derive a probable distance range of 
3.4-3.6 kpc. 

It is reassuring that the distance limits derived from 
the long-term flux history are consistent with our other 
measurements, despite the variations in m. This is par- 
ticularly relevant to the transient accretion-powered pul- 
sars in which no bursts are observed. That the long- 
term averaged m varies at all in SAX J1808.4— 3658 is 
unexpected, since the timescale for the variation is much 
shorter than the characteristic evolution time of these 
systems (of order 10^ yr). The similarity of the lower 
distance limit derived from the long-term flux and the 
estimated upper limits from the bursts suggests that the 
mass-transfer rate has already reached « Mgr, and thus 
could be expected to remain at least at that level for suc- 
cessive outbursts. Alternatively, the observed variation 
may be part of a longer-term cycle, and m may recover 
in the future. 

As a caveat regarding the burst fitting procedure, we 
point out that we neglect possible systematic contribu- 
tions, which are difficult to quantify. The most dramatic 
potential effect would be if we have missed every second 
burst in the regular Earth occultations, so that the actual 
burst rate is twice what we infer. However, we can rule 
this out fairly confidently, since the required m would be 
a factor of ~ 2 higher, in which case (unless the persistent 
flux is preferentially beamed away from the line of sight) 
the distance to the source would be around 40% greater. 
This would be significantly outside the probable range 
derived here or any r an ge quoted in the literature (e.g. 
lin 't Zand et all 1200 Ij) . iWang et alJ l|200H) derived dis- 
tance contours based on the measured Ay of the optical 
counterpart, and estimated an inclination of cos« — 0.8 
for a distance of 3 kpc. For our estimated lower distance 
limit of 3.4 kpc, we infer cosi > 0.8, so that i < 30°. 

Our best-fitting models give a value Qt ~ 0.3 MeV 
per nucleon for the flux from deep in the star that 
heats the accumulating fuel layer. For an accretion rate 
« 0.06 mEdd (eq. EI), the corresponding luminosity is 
AnR'^mQb « 2 x 10^* ergs s~^. This value is three or- 
ders of magnitude la rger than the quiescent luminosity of 
SAX J1808.4-3658. iCampana etaP pOOl modeled the 
quiescent spectrum as a blackbody plus power law com- 
ponent, and obtained an upper limit to the blackbody 
flux of 2 X 10~^^ ergs cm~^ s~^. For a distance of 3.5 kpc, 
the implied X-ray luminosity is < 3 x 10^° ergs s~^. For 
a spectral model which included a hydrogen atmosphere 
model rather than blackbody, the upper limit was a fac- 
tor of two lower. This large discrepancy between the 
outwards flux deduced from X-ray burst properties and 
the flux measured in quiescence is similar to the case of 
KS 1731—260. In that case, the observed quiescent lu- 
minosity is « 3 x 10^"^ ergs s~^ (for an assumed distance 
of 8 kpc; Rutledge et al. 2002). However, KS 1731-260 
showed a super burst during outburst, for which a lumi- 
nosity 20 times greater, > 6 x lO^** ergs s~^, is required in 
order to achieve unstable carbon ignition at the correct 
depth (Figure 20 of Gumming et al. 2006). The meaning 
of these discrepancies is not yet clear, but may indicate 
some additional heating of the deep ocean that is not 
included in current models of accreting neutron stars. 
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Our new distance estimate for SAX J1808.4— 3658 
eases the constraints on ne utron star interio r mode ls from 
the quiescent luminosity. iGamoana et al.1 l)2002|) noted 
that the low quiescent luminosity implies that the neu- 
tron star core is cold, presumably indicating a massive 
neutron star core with fast neutrino emission. Yakovlev 
et al. (2004) further found that only models with hyper- 
ons in the core could explain the low luminosity. How- 
ever, they assumed a distance of 2.5 kpc, whereas our 
new distance of « 3.5 kpc implies a luminosity two times 
greater than that assumed by Yakovlev et al. (2004). 
This is enough to relax the constraints so that core mod- 
els without hyperons are viable (see Figure 6 of Yakovlev 
et al. 2004). 

The method for comparing the observed and predicted 
burst properties may also be applicable to other bursters, 
although if the accretion rate exceeds a few percent 
TTiEdd it would be necessary to adopt a more detailed 
ignition model ta king into account the rmal inertia ef- 
fects, like that of iWooslev et alJ l|2004l) . We explicitly 
assume that the proportionality between the X-ray flux 
and the accretion rate does not change with time. This 
is likely true for the part of the 2002 October outburst 
of SAX J1808.4— 3658 covered by our simulations, given 
the monotonic relationship between the X-ray intensity 
and the c haracteristic quasi-period ic oscillation frequen- 
cies (e.g. Ivan Straaten et ani2005D : but is not in gen- 
eral true for other LMXBs. It is also important to ver- 
ify the assumptions, particularly complete fuel consump- 
tion; such conditions are likely the case for very regular 



bursts such as those typically observed in GS 1826—24 
ijGallowav et alJl2004D . but are almost certainly not the 
case generally. We note that the degeneracy for the ac- 
ceptable solutions between Xq and Zq^q arises from the 
similarity between the H-fraction at ignition for the four 
bursts, which in turn results from the relatively narrow 
range of m over which the bursts in SAX J1808.4— 3658 
are observed. If burst observations over a wider range 
of m are available, we may very well be able to break 
this degeneracy and measure the full set of compositional 
properties of the burst source. 

To our knowledge, this is the first time that a se- 
quence of bursts in the pure-He ignition regime has been 
confirmed. At first inspection, the lightcurves of the 
burst s agree well w ith the time-dependent models of 
Woo slev et al.l l|2004l) , it will be interesting in future work 
to compare the observed and theoretical lightcurves in 
more detail. 
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